Classical molecular dynamics simulation identifies catechingallate as a promising antiviral polyphenol against MPOX palmitoylated surface protein.

Cumulative global prevalence of the emergent monkeypox (MPX) infection in the non-endemic countries has been professed as a global public health predicament. Lack of effective MPX-specific treatments sets the baseline for designing the current study. This research work uncovers the effective use of known antiviral polyphenols against MPX viral infection, and recognises their mode of interaction with the target F13 protein, that plays crucial role in formation of enveloped virions. Herein, we have employed state-of-the-art machine learning based AlphaFold2 to predict the three-dimensional structure of F13 followed by molecular docking and all-atoms molecular dynamics (MD) simulations to investigate the differential mode of F13-polyphenol interactions. Our extensive computational approach identifies six potent polyphenols Rutin, Epicatechingallate, Catechingallate, Quercitrin, Isoquecitrin and Hyperoside exhibiting higher binding affinity towards F13, buried inside a positively charged binding groove. Intermolecular contact analysis of the docked and MD simulated complexes divulges three important residues Asp134, Ser137 and Ser321 that are observed to be involved in ligand binding through hydrogen bonds. Our findings suggest that ligand binding induces minor conformational changes in F13 to affect the conformation of the binding site. Concomitantly, essential dynamics of the six-MD simulated complexes reveals Catechin gallate, a known antiviral agent as a promising polyphenol targeting F13 protein, dominated with a dense network of hydrophobic contacts. However, assessment of biological activities of these polyphenols need to be confirmed through in vitro and in vivo assays, which may pave the way for development of new novel antiviral drugs.

Spectroscopic insights with molecular docking and molecular dynamic simulation studies of anticancer drug 5-Fluorouracil targeting human pyruvate kinase m2.

This study was conducted to test the efficacy of 5-fluorouracil (5-FU) as an anticancer drug against the human pyruvate kinase isozyme M2 (PKM2) using spectroscopic, molecular docking and molecular dynamic simulation studies. PKM2 fluorescence quenching studies in the presence of 5-FU performed at three different temperatures indicates dynamic quenching processes with single-set of binding (n ≈ 1) profile. The biomolecular quenching constants (kq) and the effective binding constants (Kb) obtained are shown to increase with temperature. The calculated enthalpy (ΔH) and entropy changes (ΔS) are estimated to be -118.06 kJ/mol and 146.14 kJ/mol/K respectively, which suggest the possible mode of interaction as electrostatic and hydrogen bonding. Further, these values were used to estimate the free energy changes (ΔG) and that increases with temperature. The negative ΔG values clearly indicates spontaneous binding process that stabilizes the complex formed between 5-FU and PKM2. Far-UV CD spectra of PKM2 in the presence of 5-FU shows decrease in α-helix contents which point towards the destabilization of secondary structure that weakens the biological activity of PKM2. The intrinsic fluorescence study and circular dichroism (CD) spectra showed minor conformational changes of PKM2 in the presence of 5-FU. Additionally, the results obtained from molecular docking and all-atom molecular dynamic simulation study supports the insight of the spectroscopic binding studies, and strengthens the dynamic stability of the complex between 5-FU and PKM2 through H-bonding. This study establishes a paradigm of 5-FU-PKM2 complexation and the efficacy of 5-FU that compromises the biological activity of the targeted PKM2.Communicated by Ramaswamy H. Sarma.

Lipid and cholesterols modulate the dynamics of SARS-CoV-2 viral ion channel ORF3a and its pathogenic variants.

SARS-CoV-2 accessory protein, ORF3a is a putative ion channel which immensely contributes to viral pathogenicity by modulating host immune responses and virus-host interactions. Relatively high expression of ORF3a in diseased individuals and implication with inflammasome activation, apoptosis and autophagy inhibition, ratifies as an effective target for developing vaccines and therapeutics. Herein, we present the elusive dynamics of ORF3a-dimeric state using all-atoms molecular dynamics (MD) simulations at µ-seconds scale in a heterogeneous lipid-mimetic system in multiple replicates. Additionally, we also explore the effect of non-synonymous pathogenic mutations on ORF3a ion channel activity and viral pathogenicity in different SARS-CoV-2 variants using various structure-based protein stability (ΔΔG) tools and computational saturation mutagenesis. Our study ascertains the role of phosphatidylcholines and cholesterol in modulating the structure of ORF3a, which perturbs the size and flexibility of the polar cavity that allows permeation of large cations. Discrete trend in ion channel pore radius and area per lipid arises the premise that presence of lipids might also affect the overall conformation of ORF3a. MD structural-ensembles, in some replicates rationalize the crucial role of TM2 in maintaining the native structure of ORF3a. We also infer that loss of structural stability primarily grounds for pathogenicity in more than half of the pathogenic variants of ORF3a. Overall, the effect of mutation on alteration of ion permeability of ORF3a, proposed in this study brings mechanistic insights into variant consequences on viral membrane proteins of SARS-CoV-2, which can be utilized for the development of novel therapeutics to treat COVID-19 and other coronavirus diseases.

Molecular docking and MD simulations reveal protease inhibitors block the catalytic residues in Prp8 intein of Aspergillus fumigatus: a potential target for antimycotics.

Resistance to azoles and amphotericin B especially in Aspergillus fumigatus is a growing concern towards the treatment of invasive fungal infection. At this critical juncture, intein splicing would be a productive, and innovative target to establish therapies against resistant strains. Intein splicing is the central event for the activation of host protein, essential for the growth and survival of various microorganisms including A. fumigatus. The splicing process is a four-step protease-like nucleophilic cascade. Thus, we hypothesise that protease inhibitors would successfully halt intein splicing and potentially restrict the growth of the aforementioned pathogen. Using Rosetta Fold and molecular dynamics simulations, we modelled Prp8 intein structure; resembling classic intein fold with horse shoe shaped splicing domain. To fully comprehend the active site of Afu Prp8 intein, C1, T62, H65, H818, N819 from intein sequences and S820, the first C-extein residue are selected. Molecular docking shows that two FDA-approved drugs, i.e. Lufotrelvir and Remdesivir triphosphate efficiently interact with Prp8 intein from the assortment of 212 protease inhibitors. MD simulation portrayed that Prp8 undergoes conformational change upon ligand binding, and inferred the molecular recognition and stability of the docked complexes. Per-residue decomposition analysis confirms the importance of F: block R802, V803, and Q807 binding pocket in intein splicing domain towards recognition of inhibitors, along with active site residues through strong hydrogen bonds and hydrophobic contacts. However, in vitro and in vivo assays are required to confirm the inhibitory action on Prp8 intein splicing; which may pave the way for the development of new antifungals for A. fumigatus.Communicated by Ramaswamy H. Sarma.

Illuminating the structural basis of human neurokinin 1 receptor (NK1R) antagonism through classical all-atoms molecular dynamics simulations.

Advances in structural biology have bestowed insights into the pleiotropic effects of neurokinin 1 receptors (NK1R) in diverse patho-physiological processes, thereby highlighting the potential therapeutic value of antagonists directed against NK1R. Herein, we investigate the mode of antagonist recognition to discern the obscure atomic facets germane for the function and molecular determinants of NK1R. To commence discernment of potent antagonists and the conformational changes in NK1R, induced upon antagonist binding, state-of-the-art classical all-atoms molecular dynamics (MD) simulations in lipid mimetic bilayers have been utilized. MD simulations of structural ensembles reveals the involvement of TM5 and TM6 in tight anchoring of antagonists through a network of interhelical hydrogen-bonds, while, the extracellular loop 2 (ECL2) governs the overall size and nature of the pocket, thereby modulating NK1R. Consistent comparison between experiments and MD simulation results discerns the predominant role of TM3, TM4, and TM6 in lipid-NK1R interaction. Correlation between hydrophobic index and helicity of TM domains elucidates their importance in maintaining the structural stability in addition to regulating NK1R antagonism. Taken together, we anticipate that our computational study marks a comprehensive structural basis of NK1R antagonism in lipid bilayers, which may facilitate designing of new therapeutics against associated diseases targeting human neurokinin receptors.

Dissecting the Molecular Basis of Host Leucine-Rich Repeat Containing 15 Mediated Interaction with Receptor Binding Domain of SARS-CoV-2 Spike Protein: A Computational Approach.

The detection of leucine-rich repeat containing 15 (LRRC15) as a connecting link with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) underscores the possibility of its involvement in differential restriction activity of SARS-CoV-2 pathways. However, the structure-function mechanism of LRRC15 involving the receptor binding domain (RBD) of the SARS-CoV-2 spike protein and their mode of interaction is largely unknown. Using state-of-the-art AlphaFold2 and all-atom molecular dynamics simulations, our findings provide evidences of alternative binding modes of RBD with LRR units of LRRC15 having varied affinities. Contribution of both the receptor binding regions in RBD, including receptor binding motif in accommodating the LRR domain, towards the C-terminal region, emphasizes its differential role in modulating host cell receptiveness for SARS-CoV-2, the innate immune system, as well as antiviral tone. However, further experimental validations are necessary for unravelling the unknown mechanism and distinctive features of this host receptor in the COVID-19 pandemic, involving both the transmembrane as well as cytoplasmic domain.

Exploiting the potential of natural polyphenols as antivirals against monkeypox envelope protein F13 using machine learning and all-atoms MD simulations.

The re-emergence of monkeypox (MPX), in the era of COVID-19 pandemic is a new global menace. Regardless of its leniency, there are chances of MPX expediting severe health deterioration. The role of envelope protein, F13 as a critical component for production of extracellular viral particles makes it a crucial drug target. Polyphenols, exhibiting antiviral properties have been acclaimed as an effective alternative to the traditional treatment methods for management of viral diseases. To facilitate the development of potent MPX specific therapeutics, herein, we have employed state-of-the-art machine learning techniques to predict a highly accurate 3-dimensional structure of F13 as well as identify binding hotspots on the protein surface. Additionally, we have effectuated high-throughput virtual screening methodology on 57 potent natural polyphenols having antiviral activities followed by all-atoms molecular dynamics (MD) simulations, to substantiate the mode of interaction of F13 protein and polyphenol complexes. The structure-based virtual screening based on Glide SP, XP and MM/GBSA scores enables the selection of six potent polyphenols having higher binding affinity towards F13. Non-bonded contact analysis, of pre- and post- MD complexes propound the critical role of Glu143, Asp134, Asn345, Ser321 and Tyr320 residues in polyphenol recognition, which is well supported by per-residue decomposition analysis. Close-observation of the structural ensembles from MD suggests that the binding groove of F13 is mostly hydrophobic in nature. Taken together, this structure-based analysis from our study provides a lead on Myricetin, and Demethoxycurcumin, which may act as potent inhibitors of F13. In conclusion, our study provides new insights into the molecular recognition and dynamics of F13-polyphenol bound states, offering new promises for development of antivirals to combat monkeypox. However, further in vitro and in vivo experiments are necessary to validate these results.

Machine learning and classical MD simulation to identify inhibitors against the P37 envelope protein of monkeypox virus.

Monkeypox virus (MPXV) outbreak is a serious public health concern that requires international attention. P37 of MPXV plays a pivotal role in DNA replication and acts as one of the promising targets for antiviral drug design. In this study, we intent to screen potential analogs of existing FDA approved drugs of MPXV against P37 using state-of-the-art machine learning and computational biophysical techniques. AlphaFold2 guided all-atoms molecular dynamics simulations optimized P37 structure is used for molecular docking and binding free energy calculations. Similar to members of Phospholipase-D family , the predicted P37 structure also adopts a ß-α-ß-α-ß sandwich fold, harbouring strongly conserved HxKxxxxD motif. The binding pocket comprises of Tyr48, Lys86, His115, Lys117, Ser130, Asn132, Trp280, Asn240, His325, Lys327 and Tyr346 forming strong hydrogen bonds and dense hydrophobic contacts with the screened analogs and is surrounded by positively charged patches. Loops connecting the two domains and C-terminal region exhibit high degree of flexibility. In some structural ensembles, the partial disorderness in the C-terminal region is presumed to be due to its low confidence score, acquired during structure prediction. Transition from loop to ß-strands (244-254 aa) in P37-Cidofovir and its analog complexes advocates the need for further investigations. MD simulations support the accuracy of the molecular docking results, indicating the potential of analogs as potent binders of P37. Taken together, our results provide preferable understanding of molecular recognition and dynamics of ligand-bound states of P37, offering opportunities for development of new antivirals against MPXV. However, the need of in vitro and in vivo assays for confirmation of these results still persists.Communicated by Ramaswamy H. Sarma.

Structural insights into SARS-CoV-2 main protease conformational plasticity.

The spread of different severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants underscores the need for insights into the structural properties of its structural and non-structural proteins. The highly conserved homo-dimeric chymotrypsin-like protease (3CL MPRO ), belonging to the class of cysteine hydrolases, plays an indispensable role in processing viral polyproteins that are involved in viral replication and transcription. Studies have successfully demonstrated the role of MPRO as an attractive drug target for designing antiviral treatments because of its importance in the viral life cycle. Herein, we report the structural dynamics of six experimentally solved structures of MPRO (i.e., 6LU7, 6M03, 6WQF, 6Y2E, 6Y84, and 7BUY including both ligand-free and ligand-bound states) at different resolutions. We have employed a structure-based balanced forcefield, CHARMM36m through state-of-the-art all-atoms molecular dynamics simulations at µ-seconds scale at room temperature (303K) and pH 7.0 to explore their structure-function relationship. The helical domain-III responsible for dimerization mostly contributes to the altered conformational states and destabilization of MPRO . A keen observation of the high degree of flexibility in the P5 binding pocket adjoining domain II-III highlights the reason for observation of conformational heterogeneity among the structural ensembles of MPRO . We also observe a differential dynamics of the catalytic pocket residues His41, Cys145, and Asp187, which may lead to catalytic impairment of the monomeric proteases. Among the highly populated conformational states of the six systems, 6LU7 and 7M03 forms the most stable and compact MPRO conformation with intact catalytic site and structural integrity. Altogether, our findings from this extensive study provides a benchmark to identify physiologically relevant structures of such promising drug targets for structure-based drug design and discovery of potent drug-like compounds having clinical potential.

An immunoinformatic approach towards development of a potent and effective multi-epitope vaccine against monkeypox virus (MPXV).

Monkeypox is a viral zoonotic disease, often transmitted to humans from animals. While the whole world is haggling with the COVID-19 pandemic, the emergence of the monkeypox virus (MPXV) arose as a new challenge to mankind. Till date, numerous cases related to the MPXV have been reported in several countries across the globe, but, its momentary distribution in the current time has left everyone in fright with increasing mortality and limited clinically approved treatments. Therefore, it is of immense importance to develop a potent and highly effective vaccine capable of inducing desired immunogenic responses against the highly contagious MPXV. Herein, using various immunoinformatic and computational biology tools, we made an attempt to develop a multi-epitope vaccine construct against the MPXV which is antigenic, non-allergen and non-toxic in nature and capable of exhibiting immunogenic behavior. The sequence of vaccine construct was designed using the proposed 4 MHC-I, 3 MHC-II and 4 B-cell epitopes linked with suitable adjuvant and linkers. The modeled structure of the vaccine construct was used to assess its interaction with the Toll-like Receptor 4 (TLR4) using ClusPro and HADDOCK. All-atoms molecular dynamics simulation of the MPXV vaccine construct-TLR4 complex followed by a high level of gene expression of the construct within the bacterial system affirmed its stability along with induction of immunogenic response within the host cell. Altogether, our immunoinformatic approach aid in the development of a stable chimeric vaccine construct against MPXV and needs further experimental validation for its immunological relevance and usefulness as a vaccine candidate.Communicated by Ramaswamy H. Sarma.